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Abstract 



The measurement of sin 2/3 is discussed within and beyond the standard 
model. In the presence of new physics, the angle (3 extracted from the global 
fit(denoted by fi^ff) and the one extracted from B — > J/ipKs (denoted by 
ftj/ip) are in general all different from the 'true' angle (3 which is the weak 
phase of CKM matrix element V t * d . Possible new physics effects on the ratio 
Rg = sin2/3j/^/sin2 ( 3j^ / is studied and parameterized in a most general 
form. It is shown that the ratio Rg may provide a useful tool in probing new 
physics. The experimental value of Rg is obtained through an update of the 
global fit of the unitarity triangle with the latest data and found to be less 
than unity at la level. The new physics effects on Rg from the models with 
minimum flavor violation (MFV) and the standard model with two-Higgs- 
doublet (S2HDM) are studied in detail. It is found that the MFV models 
seem to give a relative large value Rg > 1. With the current data, this may 
indicate that this kind of new physics may be disfavored and alternative new 
physics with additional phases appears more relevant. As an illustration for 
models with additional phase beyond CKM phase, the S2HDM effects on Rg 
are studied and found to be easily coincide with the data due to the flavor 
changing neutral Higgs interaction. 

PACS numbers: ll.30.-j, 12.60.-i, 13.25.-k 
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I. INTRODUCTION 



Recently, the BaBar and Belle Collaborations have reported their new results on the 
measurements of time dependent CP asymmetry Acp of decay mode B — ► J/ipKs- The 
definition of Acp is given as follows 

Acp{t) _ r - J /^ K s) - r(s°(t) - J/^K S ) 



r(B°(t) -> J/v^s) + r(5°(t) -> 

= -sin2/?j/^sin(Am jB t). (1) 

In the framework of the Standard Model (SM), the angle /3j/^ is expected to be equal to 
the angle (3 which concerns the weak phase of Cabbibo-Kabayashi-Maskawa (CKM) matrix 
element V* d 

which is also one of the angles of the unitarity triangles containing d and b quarks. The new 
measurements give the following values [0,0: 

sin2/3j/^| exp = 0.59 ± 0.14 ± 0.05 (BaBar) 

sin2/3 w | e ^ = 0.58^+of (Belle). (3) 

Combining with the ones measured earlier by CDF and ALEPH Collaboration @, 
and 0.84^04 ± 0.16 H, one arrives at the new world average value of sin 2(3 



sin 



2A7/^U = 0.59 ±0.14. (4) 



Although this result is consistent with the prediction of Standard Model (SM), the results 
from BaBar and Belle may imply a possibility that the angle measured directly from 
the time dependent CP asymmetry may be different than the one from the global fit in SM ( 
denoted by (3^) which contains the indirect CP violations of both neutral K and B meson 
mixings. This possibility has aroused many discussions in the literatures PHlO|, ^ , Hl ■ As 



has been pointed out in Ref. J8|, a lower value of sin 2(3 less than 0.5 may imply that some 
hadronic parameters are out of the reasonable range. It may imply a large bag parameter 
Bpc for the matrix element < K°\Hw\K° >, a large 577(3) breaking factor between B° and 
B° s or a small value of \V u b\. In the SM, a lower bound of sin 2(3 can be derived from the 
evaluation of ex and Am#, a conservative bound was found to be sin 2(3 > 0.42 || which is 
compared with 0.34 in |HJ due to the use of different values of hadronic parameters. The 



existence of a lower bound on sin 2(3 also holds for new physics models with minimal flavor 
violation (MFV), namely a class of models which has no new flavor changing operators 
beyond those in the SM and no additional weak phases beyond the CKM phase. 

In the SM, the angle (3^ extracted from the global fit should be the same as the one 
measured from the time dependent CP asymmetry in the decay B — ► J/ipKs- However, if 
there exists new physics beyond the SM, the situation may be quite different, the angle 
extracted from two different approaches are in general not equal. This is because ex and 
Amp, as well as B — > J/tpKs will receive contributions from new physics in a quite different 
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manner comparing to the ones from the SM. As a consequence, the extracted "/?" from two 
ways will become different. In the most general case, one may find that f3f?f ^ (3j/^, ^ (3 
when new physics exists and the ratio 

_ sin2/3j/^ 

^ = siu2^f (5) 

is therefore not equal to unity. It will be shown in detail below that the true angle (3 may 
not be directly obtained from the measurement of /3f-j? or (3j/^. 

It is obvious that the deviation from Rp = 1 provides a clean signal of new physics. 
While the different new physics models may modify its value in different ways. In the 
definition of Rp, the value of sin2(3j/^ is model dependent, the prediction of sin2/3j/</, may 
vary largely from different new physics models. On the other hand, even when sin 2/3*$? is 
extracted via the same formulae as the ones in the SM and its value only depends on the 
hadronic parameters and the experimental data, namely it looks like model independent, 
but possible existence of new physics implies that (3^$? may not be necessarily equal to the 
true (3 in the SM. This is because the experimental data should contain all contributions 
from both SM and beyond the SM if new physics truly exists. Thus the relation between 
/3f$? and (3 depends on models. In this paper, we present a general parameterization for 
possible new physics contributions to Rp. For a detailed consideration, we investigate two 
interesting and typical models, one is the model with MFV and another is the simplest 
extension of SM with just adding one Higgs doublet. For convenience of mention in our 
following discussions, we may call such a minimal extension of the standard model with two 
Higgs doublet as an S2HDM jO-IITJ, in which CP violation could solely originate from the 



Higgs potential f]T8| , pL5|,pL6|| . The experimental value of Rp is obtained from an update of the 
global fit for the unitarity triangle (UT). It will be seen that the ratio Rp may provide a 
useful tool for probing new physics. The current data lead to a low value of Rp which is not 
likely to be accommodated by the models with MFV. If the future experiments confirm such 
a low value of Rp, the models involving new interactions with additional phases beyond the 
CKM phase will be preferred. Taking the S2HDM as an example, we illustrate that the new 
physics models with additional CP phases can easily explain the current data. 

This paper is organized as follows: in section H, the ratio Rp with possible new physics 
effects is introduced and the new physics effects on Rp are parameterized in a most general 
form. In section M, the profile of the unitarity triangle is updated with the latest data and 
the resulting value of Rp is obtained. The new physics effects of the models with MFV and 
the S2HDM are discussed in section IV. Our conclusions are made in the last section. 



II. NEUTRAL MESON MIXING WITHIN AND BEYOND THE SM 

In the study of quark mixing and CP violation, it is convenient to use the Wolfenstein 
parameterization |20| in which the CKM matrix elements can be parameterized by four 
parameters A, A, p and rj. The unitarity relation 

VudV* b + v cd v; b + v td v* = o (6) 

can be depicted as a closed triangle in the p, r\ plane (see Fig.|l|). The two sides of the 
triangle can be written as: 
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R u EE = ^ = VP*+?^ 



with p ee (1 — 4^)p, and 77 = (1 — ^-)r). The three angles of a, /3 and 7 of the triangle are 
defined as follows 



sin 2a 
sin 2(3 



2r] (rj 2 + p 2 - p) 
(p 2 + 7f)((i-p) 2 + rf) 
277(1 -P) 

(i-py + v 



=2' 



sin2 7 = ^^, (8) 

p^ + T) 2 

with a + /3 + 7 = 180°. One of the important goals in the modern particle physics is to 
precisely determine those angles and to test the unitarity relation of Eq.(||). Many efforts 
have been done for that purpose and the current data have constrained the apex of the 
triangle in a small region. For example, the values of A, \V c b\ and \V u i\ are extracted from 
the semileptonic K and B decays. Besides the semileptonic decays, the important constraint 
may come from the neutral meson mixings, such as K° — K° and B? s n — B^ mixings. 

There are two important observables in neutral meson mixings: the indirect CP violating 
parameter ex of kaon and the mass difference of neutral B meson Am B(s) . In the SM, their 
expressions read 



C SM 
Z K 



{v [(1 - p)A 2 V2 S (x t ) + P e ] A 2 B K C e } e*i 



Am *t) = ^2 m w m B M (f B{3) jB^\ B S,{x t )\\^\ 2 . (9) 
where So(xi) with (xj = mf/My/) are the integral function arising from the box diagram [22 



and rjiS are the QCD corrections with the values |23| rj\ = 1.38±0.20, 772 = 0.57±0.01, 773 = 
0.47 ± 0.04 and r] B = 0.55 ± 0.01, G F = 1.16 x 10~ 5 GeV~ 2 is the Fermi coupling constant, 
f B and Bk , B are the decay constant and bag parameter for kaon and B meson respectively. 
The two constants C t and P e have the values |]S|,|l0fl: 

C e = G 2 F F 2 K m K m 2 w \ w /(<oV2ti 2 Am K ) = 0.01 

and 

Pe = [mS (x c ,x t ) - ViS (x c )] /A 4 = 0.31 ± 0.05 

In the global SM fit of UT, one usually assumes that the SM is fully responsible for 
the observed experimental data. This scenario, however may be invalid if there are new 
contributions beyond the SM. If new physics effects exist, then e e ^ p ' 7^ e^- M (p, 77) and Am^ 7^ 
Amg M (p,r]) . In general, one may write 

*T = 4 M (M)[l + 5e K ] ee e s K M {p™,f™) 
Am7 = Am s B M (p, fj) [1 + 6m B ] ee Am s B M (p s f ^, fjfj<) [1 + 5m B ] (10) 
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where Sex and 5m B ( Stub) represent possible new physics contributions. Therefore the 
parameters p and rj extracted from the fit are in general the effective ones which may 
significantly deviate from their true values in the SM . The effective values p^ and rj^ff 
are defined as follows: 



= {fflf [(1 - P^)A 2 m So(x t ) + P] A*B K C t } e *i 
Amr/[1 + Sfh B ] = g_(l - \fA^mlm B (f B VB)\ B S (x t ) ((1 - p™f + (fg ) 2 ) . (11) 



where, e e ^ p and Am^ are the values measured from the experiments | 31| . It is seen that 
in the existence of unknown new physics, one may even not be able to extract the effective 
parameters p^ 1 and fj^^f from the data e e ^ p and Am e ^ p . Only when the new physics 
contributions to the term 5m B are neglected, the two effective parameters p s ^ and fjf^f 
can be determined from the data e e £ p and Am e ^ p . In our following fit, due to the large 
uncertainties in both theoretical parameters and the present data, we actually take 5m B = 
as an approximation, so that the angle can be extracted via the same form as the one 
in Eq.(|) 



sin2/? SM = Vfit { Pfit ' (12) 



From Eq.flT2]) , by eliminating So(x t ) in the expression for ex and Am B , the following relation 
holds 

i2^' = C K (l-4 i f). (13) 

C K = t^gf^ 1 and u = C,A*B K P M . 
Before proceeding, we would like to emphasize that even if taking the expressions satis- 
fying the same relations as the ones in the SM, it is clear that the resulting values for p^ 
and rjf-f are not necessary to be the true ones in the SM if new physics beyond the SM 
truly exists. In general, one has pfjjf 7^ p and m^f 7^ 77. Only when there are no new physics 
contributions to the K° — K° and B° — B° mixings, i.e. e e ^ p = ef- M and Am e B p = Am'™ , 
one can then get the true values p s ^ = p and fjftf = fj. Thus the fit can really probe 
the true value of sin 2(3. However, if there exists new physics beyond the SM, the relation 
between sin 2(3^ and sin 2/3 will become complicated and depend strongly on new physics 
models, which can explicitly be seen below. Note that the values of pj^f and r)f*f fitted in 
such a way remain model independent though they may not be equal to the true values of 
p and r\ in the SM. 

It is noticed that in the SM neutral meson mixings only occur at loop level, new physics 
contributions are expected to be considerable. If we denote the transition matrix element 
< f°\H e ff\f° > of neutral meson /°(/° = K°, B9^) by M{ 2 , then the new physics effects on 

can be parameterized in the following general form: 

M L = E AfAjMf 2 ^(l+4) + M/ 2 (14) 

i,j=c,t 

where = c,t) are the products of CKM matrix elements with the definitions \f = 
V* d V is , Af = V* d V ib , and \f s = V*V ib . The matrix elements M™/ are the ones in SM. 
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They receive contributions from the internal tt, cc and ct quark loops in the box diagrams 
which are responsible for the f° — f° mixing. Real parameters rL reflect the possible new 
physics contributions to the corresponding loop diagrams, which carry no additional weak 
phases besides Models with such property include some versions of SUSY models and 
type I and II 2HDM. Obviously those models belong to the ones with MFV since they carry 
no additional phases. New physics contributions with new weak phases are all absorbed into 
the last term M[ 2 . This term may arise from additional new interactions, such as the tree 
level Flavor Changing Neutral Current (FCNC) and the superweak-type interactions [19] as 
well as interacting terms with new weak phases. 

In evaluating ck in kaon system, although the contribution of tt quark loop is dominant, 
the effects of cc and ct quark loops are not negligible. While in the calculations of the mass 
differences of B^ meson, one may only consider the tt quark loop in the box diagram. 

By using the Wolfenstein parameterization for CKM matrix elements and the parame- 
terization of Eq.(H), in the presence of new physics, the expression of ex is 



[n [(i - p)A%s Q (x t ) (1 + 



K 



+ PJ1+T, 



ctj 



A 2 B K C e ) 



+ £k 



(15) 



with 



7] 3 S (x c , x t )r% - rjiS (x c )r, 



„K 
cc 



T] 3 S (x c ,X t ) - T]iS (x c 



V2A 



-ImM 



m K 



K 
12 



(16) 



where the new physics effects with MFV have been expressed as the corrections to the 
integral function So(x t ) and P e . The effects with additional phases have been absorbed into 
Ik- The B^ — B?,^ mass difference is also modified as follows 



Arris 
Am Ba 



M- 
M 



Am 



SM 
B 



Am s B M 



1+rg 



e -2i P+ ~B e ^ 

l + r B "e 2i<t>I 



(17) 
(18) 



Where we have kept only the tt quark loop in the box diagrams, the contributions from 
other loops with internal cc and ct quarks are highly suppressed and can be safely neglected. 

From the definition of Acp, the (3j/^ measured from the time dependent CP asymmetry 
will be the total phase of Mj|. Thus in the presence of new physics, The angle j3j/^ extracted 
from the time dependent CP asymmetry Acp{t) in decay B — > J/ipKs is in general different 
from the true value (3 in SM. This can be seen explicitly from the following relation 



sin 2/3jyty 



sin 2(3 — r sin 20 



l + 2r B cos2(0 B + /?) + (f s ) 2 



(19) 



The usual fit of the profile of the UT has also been made under the assumption that there 
is no new physics beyond the SM. The angle u /3f-f" which is obtained from the measurement 
of K° — K° and B° — B° mixing is usually assumed to be the true angle (3 associated with 
the phase of V* d . In general, one should not make any assumption for purpose of new physics 
probing, since it may be greatly modified when new physics effects come into the evaluation 
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of 6k and Am^,,. 

the ones in SM to extract the angle (3, we may only get the effective values of pj-f , r/j-i , 



As we have discussed before even if we still use the same formulae as 

SM „SM 



and fif^ 1 . Their values may be far away from the true ones, this is because the experiment 
measured results for ex and Am# (s) are the total ones which actually include all possible 
new physics contributions if new physics truly exists. 

Combining Eq.(^) and ( O ) and using the definition of sin 2/3 in Eq.(^), we come to the 
following simple relation 



sin 2/3 = C K 



1 



,B 



K 



1 - 



— u>7)(l + r, 



ct 



2f B cos2(0 B + /3) + {f B ) 2 



(20) 



The relation between /3«^ an d P * s straight forward, 



l-ufjfg) sin 2(5 



sin 2(5 



SM 
fit 



K N 



B 



1 + 2f B cos2(0 B +/3) + (f B ) 2 [|1 - e K /e K \ - ufj(l + r ct )} V 1 + r u 



(21) 



To probe new physics effects, let us define the ratio Rp between the two observables 
sin 2/3//,/, and sin 2/3/^ as follows 



Rr 



sin 2(3 j /rp 
sin2/3f M 



1 + 



fit 

B 



K 



l + rf t 



[l - <W e /d - + r ct) 

1 ^Vfit 



, B sin 20 g 
r sin 2/3 



(22) 



Note that the factor yl + 2f B cos2(0 B + /3) + (f B ) 2 in Eq.(PD|) and (|TTf ) cancel each other 
in the expression of i?^. In the SM, Rp is clearly equal to unity. From the above equation 
the deviation from unity due to the new physics effects can be expressed as the products 
of three terms: (1) The ratio of new physics corrections to the tt quark loop in the box 
diagram between B° — B° and K° — K° mixings. Since almost all the new physics models 
give similar contribution to B° — B° and K° — K° mixing, their effect may cancel in some 
extent. (2) The corrections to the cc and ct-quark loops in the box diagram in kaon meson 
as well as the ones from the additional contributions to K° — K° mixing which contain 
new CP-violating phases. The corrections from cc, ct loops are suppressed due to the small 
quark mass m c if the additional contributions come from the new scalar particles such as 
the charged Higgs particle H ± in some models. However, in general case, the corrections 
may be considerable. Although the small mass difference between K° and K° has imposed a 
strong constraint on the real part of M\2 , its imaginary part can still be large and result in a 
sizable Ik such as the sup erweak- type interactions. (3) The corrections from the additional 
contributions to B° — B° with new CP-violating phases. It can be seen that the corrections 
depend on the product of r B and sin0 s . If <p B = the correction to Rp disappears. This is 

directly related to the cancelation of the factor \Jl + 2r B cos2(0 B + (5) + (f B ) 2 in the ratio 
Rp. Although the new physics with zero 4>b can affect both (5j/^ and (5f it , their net effect 
on Rp is diminished. 
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III. Rp AND (3 s F f T FROM THE GLOBAL FIT 

To get the value of Rp from the data, one needs to know the values of (3j/^ and Pfu- 
Unlike the measurement of (3j/^ which is theoretically clean, the extraction of (3f-jf seems to 
surfer from some uncertainties as it involves parameters with large theoretical uncertainties. 
There is a long time debate on the treatment of the theoretical errors. Basicly there are 
three different approaches [Q. 

(1) Assuming all the errors, both experimental and theoretical ones are the random 
variables which obey Gaussian distribution [p6[l . The most recent fitted results at 95%CL 
are ||25||: 



- 0.82 < sin2ag f < 0.42, 0.49 < sin2/^f < 0.94, 42° < 7 g f < 83° (23) 

(2) The distribution of experimental error is treated as Gaussian, but one imposes a 
prior flat distribution on the theoretical error, i.e. assuming the theoretical parameters are 
the random variables which are uniformly populated in some reasonable regions , then by 
using the Bayesian approach, the distributions of the parameters can be determined p4],|27|. 



The recent fits using this method give the results [24 



- 0.88 < fan2ajg < 0.04, 0.57 < sin2/^f < 0.83, 42° < jf*f < 67° (24) 

(3) The distribution of experimental error is treated as Gaussian, but one does not 
assume any distribution on the theoretical error, since they are unknown parameters which 
should have a unique value, their errors should be largely deduced with the improvement 
of the theory and the experiments. To get a quantitative result, the space of the allowed 
region for the theoretical parameters are scanned. For each set of parameters a contour at 
some confidence level ( for example 68% or 95%) is made, then the whole region enveloping 
all the contours is considered as the allowed region at some overall confidence level. Note 
that in doing this , the final results will have no clear probability explanation. The most 
recent fit gives [p8| , [2"^ 

- 0.95 < sin2ag f < 0.5, 0.5 < sin2/^f < 0.85, 40° < 7 £f < 84° (25) 

In this paper, to get a more conservative conclusion, we adopt the last method and 
update the fit with the latest data on Amg and Amg s 

Let us briefly describe the method and parameters used in the fitting of the unitarity 



triangle (UT). The detailed description of the fitting procedure can be found in Ref. POfl . 

Among the four Wolfenstein parameters A,p,rj and A, the value of A has been well 
determined with a relative high precision through semileptonic kaon decays, K + — >• 7i°e + u e 
and K\ — > 7r~e + z/ e . In the fit, we quote the value of A = 0.2196 ]3TJ] and take it as a fixed 
parameter. 

The constraints on the apex of UT can be obtained from various experiments of the 
semileptonic B decays. In b — > c transition such as B — > XJu and B° — > D* + l + i7 e , the value 
of \V c b\ can be determined. Here we quote the value of V c b from the most recent LEP results 



32] 



\Vcb\ = (40.4 ± 1.8) x 10" 3 , (A = 0.850 ± 0.037) (26) 
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which is an average of measurements between inclusive and exclusive B decays. From 
B — > pl + v e the value of \V u b\ could be extracted. However, in the determination of K 



ub • 



some model dependences have to be involved in the evaluation of the ratio [Kb/Kb I, which 
results in a considerable theoretical error. For a recent review of the theoretical error in 



determining | Kb/Kb | , we refer to the reference [35]. The recent measurements of Kb from 
LEP and CLEO collaboration give the following values: p6L|3^ 



IK 



ub 



(4.13^(stat.+det.)^ 3 48 (6 - c syst.)t^(6 
±0.02(r 6 ) ± 0.20(Model) x 10~ 3 .(LEP) 

N+0.21 



u sys. 



V ub \ = (3.25±0.14(stat.)l^(syst.) ±0.55(model) x lO^.(CLEO) 



(27) 
(28) 



Note that both the determination of \V u b\ and |Kb| may still suffer from sizable theoretical 
uncertainties pS|-p][|. In the fit, we take the central value of iKb/Kbl as a f ree theoretical 
parameter which lies in the range 



< V ub /V cb >G [0.07,0.1] 
and the value of [Kb I has the following form 

I Kb I =< Kfc/K* > |K c6 | ± 0.14 



(29) 



(30) 



It is noted that possible new physics contributions to the semileptonic bottom meson decays 
are in general expected to be small, the extracted CKM parameters |K&| and |Kb| are 
regarded as the true ones |Kb| = A\ 2 and |K&| = V P 2 + V 2 - 

The other important constraint comes from the measurements of ex- The expression of 
€k has been discussed in the previous section. From Eq. ([l5|) , the major theoretical error 
comes from the hadronic parameter Bk, and the QCD correction rji, and r] 3 . the recent 



lattice calculations give the result j|l|] 

0.8 <B K < 1.1. 



(31) 



The errors of r\i can be found in the previous section. There are also experimental errors on 
the values of c and t running quark masses, 



m c (m c ) = 1.25 ± 0.1 GeV, rh t {m t ) = 165.0 ± 5.0 GeV 



(32) 



The constraint of Rt could arises from both B° — B° and B° s — B° s mixing. Here the 



parameter with a large uncertainty is /bVBb- In the fit, we take the value [H2 



160 MeV < f B jB B < 240 MeV 



(33) 



Recently, the BaBaR and Belle Collaboration has reported their first measurements on 
ArriB d - The combined result from the hadronic and semileptonic decays is [0,0] : 



Am B d 

Am Bd 



0.512 ± 0.017(stat.) ± 0.022(syst.) ps' 1 (BaBaR) 
0.456 ± 0.008(stat.) ± 0.030(syst.) ps' 1 (Belle) 



(34) 



comparing with the old world average value of 0.472 ± 0.017ps 1 [[Hp, the result of BaBaR 
is slightly higher. Now the new world average is 
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Am Bd = 0.478 ±0.013 ps~ l (35) 

Besides Am B , the mass difference of B® and B° s also imposes a strong constraint on the 
size of R t . It is helpful to introduce a £77(3) breaking factor £ 2 = (f Bs \/B Ba ) 2 / (fB\/B B ) 2 , 
the constraint on R t from Am B reads 



* < f iF^T^f- (36) 
A V m Bd Am Ba \ min 

In the fit , we take £ G [1.08,1.2] [[44]]. The lower bound of Am Ba is now updated to 
Am Bg > 14.9 ps _1 at 95%CL []34| which is also higher than the previous value of 14.3 ps _1 . 
This bound is obtained through the so called 'amplitude method' |53). The latest data 



indicate that Am Ba ~ 17.7 ps . If this result can be confirmed by the future experiment, 
Am Bg will impose strong constraint on new physics models |[45|| . 

The other two important CP-violation observables are the direct CP violation parameter 
e'/e and sin 2(5 from the time dependent CP asymmetry of decay B — > J/ipKg. As the 
measurements of sin 2(5 and e' / e are now preliminary, they are not considered in our present 
fitting . However, it needs to emphasize that including the constraints from e'/e may impose 
a nontrivial lower bound of r\ > 0.32 [pgj] . It would be interesting to further discuss its 
constraint with the improved prediction for e'/e . 



The basic idea of the fit is based on the least square method. The first step of the 
procedure is to construct a quantity x 2 which has the following form 

^ = E (/,(A M )-</,» 2 (37) 

i a i 

where < fi > and Oi are the central values and corresponding errors of the experimentally 
measured observables. They contain | V c b\, \V u b\, \ck\, Am B and m c ,m t . fi(A,p,r]) are the 
values calculated from the theories. They are the functions of parameters A, p, 77. The set 
of A, p, 77 which minimizes the x 2 wm De regarded as the best estimated values. On the 



determination of Am Bs , the ' amplitude method' J33| is adopted, which uses the amplitude 
curves measured from the experiments 133] and add a term \ 2 S to the total x 2 : 



2 (l-A(Am Bs )f 



X s = ~ (38) 

Note that there is an alternative approach to build the log-likelihood function in which 

one uses the values relative to the ones at Am s = 00, i.e. use [ i, 1 ' 4r] instead of 

3^,2i,5|]. The detailed comparison between these two approaches can be found in Ref. [24 



In the fit, we choose A, p, r\ and m c , mt as free parameters to be fitted and scan the allowed 
range for the theoretical parameters < V u b/V c b >, Bk, f B y/B B , i] s , 771, and ^3. For each set of 
theoretical parameters, a 68%(95%) CL contour which corresponds to x 2 — Xmin + 2.4(6.0) 
is made, then the region enveloping all the contours is the allowed region for p and rj at an 
overall 68% (95%) level. As usual, a x 2 probability cut is used to reject the contours with 
relative high Xmin which means that the fit is not consistent. 



The fit is implemented by using the program package MINUIT 47 . The result is shown 



in Fig.|[ From the figure, one can read off the allowed region for p and 77. The results are 
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0.09 < p% [ < 0.28 , 0.23 < rjffl < 0.46 (at 68% CL) 

0.10 < p s f fl < 0.30 , 0.20 < rjfg < 0.49 (at 95% CL) (39) 

or equivalent ly, 

-0.89 < sin 2a s f */ < 0.37, 0.56 < sin 2/3^f < 0.85, 41° < sin 2 7^ < 75°, (at 68% CL) 
-0.95 < sin 2ag f < 0.41, 0.48 < sin 2/?Jjf < 0.88, 38° < sin 2 7^ < 76°, (at 95% CL) (40) 

and the range for the combination factor of the CKM matrix elements, Im\ t A 2 X 5 77, used 
in the calculation of direct CP- violating parameter e'/e in kaon decay [}46~ |, is given by 



0.76 x 10~ 4 < Im\ t = A 2 \ 5 rj < 1.73 x 10~ 4 , (at 68% CL) 

0.66 x 10~ 4 < Im\ t = A 2 X 5 r] < 1.84 x 10" 4 , (at 95% CL) (41) 

As the allowed ranges for theoretical errors which are the main sources of total errors 
are the same in both fits , the results at 68% and 95%CL are quite similar, It is needed to 
note that by using the scanning method, there is no probability explanation of the contours. 
One can not get the usual central value of the fitted parameters. 

In the fitting, the x 2 probability cut is set to be Prob(x 2 ) < 5%. The final results depend 
on the cut. We have checked that a lower cut such as 2% or 1% will give a larger allowed 
range, which allows the angle 7 to be greater than 90°. The probability of large 7 has 
been widely discussed in Refs.. |48| to meet the recent CLEO data on hadronic charmless 
B decays. However, the model independent analysis show that in general there are two 
solutions for 7 from the charmless B decays. The one with 7 < 90° and the other one with 



7 > 90° [i9|. The data of Am^ especially on Am^ may strongly constrain 7 to be less 
than 90°. From Fig.|3], with the cut of 5%, there is no indication of 7 > 90°. 

With the angle /3fu obtained from the above fit and the average value for the ratio 
R/3 is found to be 

0.57 < Rf p - < 1.1 (at 68% CL) , 0.31 < Rf p - < 1.36 (at 95% CL) (42) 

The present data prefer a small value of R@ < 1 at 1-a level. If it is confirmed by the future 
experiments, it will of course be a signal of new physics. Further more, Rp may be used to 
distinguish different effects of new physics. This will be discussed in the next section. 



IV. NEW PHYSICS EFFECTS ON THE RATIO R p 

The contributions of new physics to Rp depend on the models. In this section, we choose 
two type of new physics models to illustrate new physics effects, one of them is the model 
with MFV, the other is the models with new CP- violating phases, such as the S2HDM. 

A. Rp in models with MFV 

In a class of the models with MFV, the expression of Rp may be greatly simplified. As 
MFV implies that there are no new operators beyond the SM, the values of rf t and rf t 
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which come from the same internal tt quark loop in B° — B° and K° — K° mixings must be 
the same f9.10.21]. MFV also means that new CP- violating phases must vanish. Thus the 



following conditions hold 

r% = r£, f B = 0, <P B = 0, and ~e K = (43) 



Thus one has 



sin 2(3 j ft = sin 2/3, (44) 



which shows that in this case the angle measured from B — > J/ipKs is the 'true' angle 
(3 up to a four fold ambiguity. Assuming that r ct is small, which is a good approximation 



for many models | 21j| , the ratio Rp has the following simple form 



R p\mfv ~ 1 —§M ( 45 ) 

Since the value of \V u b\ is extracted from the semileptonic B — > n(p)ev e decay dominated 
by the tree diagrams, it is widely believed that its value is not likely to be affected by new 
physics. This may result in the relation: p 2 + fj 2 = (pf^f) 2 + {V/u) 2 = Combining 
this relation with Eqs.(|44j) and (§), the values of p and fj can be solved as functions of two 
obserbales R u and (3jft. 

p = R u cos <p±, fj = R u sin ip±. (46) 

with ip± = sin _1 (sin/3j/^/_R u ) ± (3jft. Note that the solution for ip + corresponds to p < 0. 
This is quite different from the value of pfjjf which seems be positive under the constraint 
from the data of B° s — B® mixing. If one assumes that the difference between true p and 
the fitted pfjf caused by new physics effects is not too large, the solution for ip + could be 
ignored. Taking the value of R u and sin 2/3//^,, one may find at la level 

1.03 < R/bImfv < 1-36 for ^_ . (47) 

For ip + solution, one gets a small value for R^\mfvj i- e - ; 0.90 < Rp\uFV < 1.21. In 
obtaining these numerical values, the errors of u and m^f are considered to be independent 
for simplicity. In general, they may be treated to be correlated as they all depend on the 
hadronic parameters. For a more detailed analysis, it should be interesting to include the 
correlation effects. 

Comparing R^mfv with the experimental values of R e ^ p ' in Eq.(fE^) , it is found that 
the models with MFV prefer to give a large value of Rp which is not likely to agree with 
the current data. If the disagreement is confirmed by the future more precise data, all the 
models with MFV will be ruled out. It implies that new physics effects with flavor changing 
interactions beyond the CKM quark mixing must be involved to explain the data. 



B. R p in S2HDM 



The S2HDM without imposing discrete symmetries |lT|-[r7|,|5l 
models which can give nonzero value of ix, r B and <p B . 



is a good example for 



In the S2HDM, the tree level flavor 
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changing neutral current can be induced from the scalar interactions between neutral Higgs 
bosons and quarks. As all the Yukawa couplings in the model are allowed to be complex, 



the values of r B , r B and <p B in S2HDM are given by 



the S2HDM can give rich sources of CP violation p~5|JT6|j5^1 . It is not difficult to find that 

tt ) 



r tt - ^ U yt— B ww {xt) + 2 I&I V* VB B ww {xt) W 
■ ? ,b -Bb^( 2V3nvm B \ 2 m d 1 Im(l^ 13 ) 2 

,m ^ "B B t\ m Hom t ) m b V t 2 dVB BWW(x t ) 1 j 

where £ q is the Yukawa coupling constant between charged Higgs and q quark. £>ww(hw,hh) 



f B sin : 



are the integral functions of the box diagrams [50 with variable x t = rrqjvn^^ y t = mj/m H 



and Dw = m^/raj]. i]^ H is the QCD correction. It is seen that the value of r B is directly 
related to the coupling £ t . In Fig.[|, r B is plotted as a function of the charged Higgs mass 
mn± with different values of ^. It can be seen that the typical value of r B is around 0.2. 
For large & , r B could become large. 

The value of f B depends on the imaginary part of K fc ^ which is given by reduced Yukawa 
couplings between the fc-th neutral Higgs boson and quarks. Bb and Bb are the bag pa- 
rameters for (S + P) 0(5' — P) and (V — A) (g)(V — A) four quark operators respectively. 

Similar calculations can be made for the K° — K° mixing. r^,r ct in S2HDM are given 

by 

rl = r B t (50) 

VM-c&iV% H Ut\ 2 B» H (y c , y t ) + 4y» w Re(Ut)B BW (y c , y t , y w ) 
Tct (%)r lB BWW(x c ,x t ) - r h B»'W(x c ) 1 ' 

(52) 

Since x t /x c is of order (9(10 4 ), can be approximately written as 



r 



« ~ I / — B$ H (y c ,y t ) u ( 2 

at ~ 2 VVtVc ^ B WW( XcjXt 



:\ Ut\ 2 (53) 



As the value of y/yty c is of order 10~ 2 and |£ c 6| 2 is typicly of order 1, one may find that r ct is 
relatively small and of order (9(10 -2 ~ 10 _1 ). The S2HDM contribution to e# is complicate, 
it may arise from both short and long distance interactions. The detailed discussion of Ik 
can be found in Ref. [|17[ 



If one ignores the r^ t , the ratio Rp mainly depends on the value of , ImY^ 13 and 
m#o as well as the true value of sin 2/3 and f/. The ratio Rp has the following form 



Rr 



[1 ~ e K /e K \ -ur] 
1 _ um SM 



' B B ^ ( 2V37ivm B \ 2 m d 1 lm(Y k d 13 f 
^ B B ^{ m H om t ) m b \V*\ VB B ww (x t ) ' U)) 



Which shows that the new interactions between neutral Higgs and quarks are the main 
sources for the changing of ratio Rp. The value of Rp as the function of neutral Higgs mass 
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m#o is plotted in Fig.[5]with different values of ImY^ l3 and ex/^K- Where for simplicity we 
have taken fj ~ fjfn = 0.35 and Bb — Bb- From the figure, the value of Rp can be smaller 
than 1 for positive Im(Y^ 13 ) 2 or ek/ek- Thus the present data can be easily explained in 
the S2HDM. 



V. CONCLUSIONS 

In conclusion, we have investigated the measurement of the UT angle f3 within and 
beyond the SM. It has been shown that if new physics beyond the SM truly exists, the angle 
P extracted from the global fit (denoted by /3f# ) and the one extracted from B — ► J/tpKg 
(denoted by could be in general all different from the 'true' angle (3 in the SM. The 

observable Rp, which is the ratio between the angle f3 measured from the time dependent CP 
asymmetry of decay B — > J/ipKg and the one extracted from the global fit, is introduced 
and studied in detail. By using the scanning method, the value of Rp is obtained from an 
update of the fit of the UT with the latest data and it is found to lie in the range. 

0.57 < Rp xp - < 1.1 {at 68% CL) , 0.31 < Rj p - < 1.36 {at 95% CL) 

As the SM calculation gives Rp = 1, the deviation of Rp from unity is a clean signal of new 
physics. On the study of new physics effects, we have given a general parameterization of 
Rp. By using this parameterization, the effects from the models with MFV and the S2HDM 
have been investigated in detail. It has been found that the MFV models seem to give a 
large Rp relative to the current data. If it can be comfirmed by the future experiments, 
it will indicate that this kind of new physics are difficult to explain the small value of Rp. 
Such models include type I and type II 2HDM and some simple versions of SUSY. Thus new 
physics with additional phases should be considered. The most recent discussions on models 
with MFV. can also be found in Ref. • As an illustration for those models with 

additional phases beyond the one in the CKM matrix element, the S2HDM effects on Rp 
has been studied, especially, if the superweak type contributions fl9| to ek are considered, 
it could easily coincide with the data. 
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FIGURES 
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FIG. 1. The unitarity triangle 
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FIG. 2. The new physics contribution to M\2 



17 
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(b) 

FIG. 3. The allowed region for p and r\ from the global fit at 68%CL (a) and 95%CL (b). The 
individual constraints from V u b, ex and B9 n — §9, are also shown. 
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FIG. 4. rf t as a function of the charge Higgs mass m H ±. The three cures correspond to 
&=0.6(solid),0.4(dash),0.2(dotted) respectively. 
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FIG. 5. Rp as a function of the mass of neutral Higgs mass mh with different values of Im(Yf 13 ) 2 
and e~k/eK- The shadowed band indicates the allowed region for Rp at 68% CL. 

(a) Jm(y 1 d 13 ) 2 =0.01(solid), 0.02 (dashed), -0.01 (dotted), -0.02 (dotted-dashed) with 
\l-e K /e K \ =0.9. 

(b) The same as (a) with \l — e~K / £r\ = 0.8. 

(c) The same as (a) with |1 — ex/ex| = 1.1- 

(d) The same as (a) with |1 — ex/ex| = 1.2. 
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